Phase transitions in a spinless, extended Falicov-Kimball model on the triangular 

lattice 
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A numerical diagonalization technique with canonical Monte-Carlo algorithm is used to study the 
phase transitions from low temperature (ordered) phase to high temperature (disordered) phase of 
spinless Falicov-Kimball model on a triangular lattice with correlated hopping (t 1 ). It is observed 
that the low temperature phases (i.e. regular, bounded, striped and segregated) persist up to a 
finite critical temperature (T c ). In addition, we observe that the critical temperature decreases with 
increasing the correlated hopping in regular, striped and bounded phases whereas it increases in the 
segregated phase. Single and multi peak patterns seen in the temperature dependence of specific heat 
(C v ) and charge susceptibility (x) for different values of parameters like on-site Coulomb correlation 
strength (U), correlated hopping (t') and filling of localized electrons (n/) are also discussed. 

PACS numbers: 



I. INTRODUCTION 

The study of correlated systems like transition-metal 
dichalcogcnidcai~— , cobaltates^, Gdl?r& and its doped 
variant Gdl^U^ have attracted considerable attention 
recently as they exhibit a number of remarkable cooper- 
ative phenomena, like valence and metal insulator tran- 
sition, charge and magnetic order, excitonic instability 1 
and possible non-Fermi liquid states^. It has been sug- 
gested recently that these correlated systems may very 
well be described by an effective Falicov-Kimball model 
(FKM) on a triangular lattice^. These systems are geo- 
metrically frustrated. The frustration gives rise to a large 
degeneracy at low temperatures and competing ground 
states close by in energy. A consequence of this is a fairly 
complex ground state magnetic phase diagram^ and the 
presence of soft local modes strongly coupled to the con- 
duction electrons^. 

The usual FKM^ was proposed to study the 
semiconductor-metal transition in SttiBq and other cor- 
related electron systems. The FKM considers only the 
kinetic energy of the itinerant electrons and the local 
Coulomb interaction between the itinerant and localized 
electrons. The effective interactions in the FKM are 
mediated by the itinerant electron a 10 ' 11 . Therefore, the 
FK model can be used to model crystallization in the 
solidsii. This model has been applied to different types 
of lattices (bipartite/non-bipartite) and in different di- 
mensions. More results are available for the FKM on 
bipartite lattice o 12-17 than on non-bipartite lattices. 

Recently several new ground states were observed for 
an extended FKM on a triangular lattice where the effect 
of correlated hopping t', was taken into account^. As we 
know FKM is applicable to the systems where both itin- 
erant as well as localized electrons are present. In the fol- 
lowing, we will denote itinerant electrons as d-elcctrons 
and localized electrons as /-electrons. The correlated 
hopping term is influenced by the number of localized 
electrons present on the neighboring sites as explained 
below. For example, in some rare earth compounds (es- 



pecially the mixed- valence compounds), the rare earth 
ions with two different ionic configurations f n and f n ~ l 
have different ionic radii. The ions with Z™" 1 configu- 
ration have less screening of nuclear charge as compared 
to the ions with f n configuration. The d-orbitals in ions 
having j n ~ x configuration feel more attraction due to the 
nucleus and hence get contracted. This leads to increased 
localization of those orbitals. Hence the d-orbital over- 
lap between nearest neighbors depends on local /- elec- 
tron occupation of neighboring ions, resulting in a corre- 
lated hopping of d-electrons. Such a correlated hopping 
term also appears in the first principles calculation^ of 
the tight binding Hamiltonian and is usually neglected 
in Hubbard type models for being smaller compared to 
the on-site term. Its significance in superconductivity 
has been emphasized alread y 20 ' 21 and in the context of 
FKM, it has been considered by several authors 2 ^— us- 
ing various approximations. 

We, therefore, extend the model to include a corre- 
lated hopping term in the Hamiltonian for the itinerant 
electrons. 

H = — + fiS^dU, +E f J2 fUi 

(ij) i 

+Uj2flfi4di + E4(/M + fUMdj- (!) 

here d\, di (fj , fi) are, respectively, the creation and anni- 
hilation operators for d-electrons (/-electrons) at the site 
i and fi is the chemical potential. The first term in Eq.(l) 
is a measure of kinetic energy of d-electrons on a triangu- 
lar lattice. The = t, if i and j are nearest neighboring 
sites and tij = 0, otherwise. The second term represents 
the dispersionless energy level Ef of /-electrons while the 
third term is the on-site Coulomb repulsion between d- 
and /-electrons. The last term (i.e. the correlated hop- 
ping term) is an extension of the usual FKM as discussed 
above. 

A comprehensive study of the ground state phase dia- 
gram of spinless Falicov- Kimball model on a triangular 
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lattice with correlated hopping at different commensu- 
rate fillings of localized electrons in (n/ — t') plane is 
reported elsewhere 1 ^. With correlated hopping (t' € 
[—1,1]) two different phases are observed at i-filling of 
localized electrons namely (i) regular and (ii) segregated 
phases for U = 1, 3 and 5. At i-filling of localized elec- 
trons (i) bounded (ii) axial striped and (iii) segregated 
phases for U = 1 and (i) striped and (ii) segregated 
phases for [7=5 are reported. Some of these charge 
ordered phases are also observed experimentally in the 
above mentioned systems. Now, question arises, what is 
the temperature stability of these phases and which type 
of the phase transition occurs from the low temperature 
ordered phase to the high temperature disordered phase. 
In this work, therefore, we study the half filled FKM in 
its spinless version for different range of interactions with 
correlated hopping f', on a triangular lattice numerically 
and observe the stability of these phases and type of the 
phase transition at finite temperature. 



II. METHODOLOGY 

In the extended Falicov-Kimball Hamiltonian (1), 
hybridization between d— and /—electrons is absent. 
Therefore, the local /-electron occupation number fif^ = 
fifi commutes with H. The Wj = fjfi is a good quan- 
tum number taking values either 1 or according to the 
site i is occupied or unoccupied by f-electron. The gauge 
symmetry also implies^ that interband excitonic aver- 
ages of the type (d\fj) are identically zero at any finite 
temperature (implying an absence of homogeneous mixed 
valence) and the /-electron level remains dispersionless 
((fifj) = 0, for i 7^ j). Following the local conservation 
of /-electron number the Hamiltonian may be written as, 

H = £ hij {u)4di + (2) 

m 

where h tj (u) = [-t tj + t'^oJi + loj)] + (Uu>i - //)<%. 

We work at half- filling, i.e., Nf + Nd = N, where 
Nf, Nd are the total number of /- and d-electrons and 
N = L 2 (L = 12) is the number of sites. For a lattice 
of N sites, the Hamiltonian H is now an N x N matrix 
for a fixed configuration u. Wc set the value of Uj=l for 
nearest neighboring sites i and j and zero, otherwise. We 
choose a particular value of Nf (0 < Nf < N), and a con- 
figuration uj = {u)x,U32, ■ ■ ■ , <-<Jn} for Nf. Choosing values 
for t' and U, the eigenvalues Aj of H(uj), are calculated 
by numerical diagonalization with the periodic boundary 
conditions (PBC). Average value of the physical quanti- 
ties is calculated by Monte Carlo sampling method. The 
details of the method is reported elsewher o 18 ' 26 " — . 

The partition function is written as 

Z = U( E Tre-W^A (3) 

i \wi=0,l / 



where the trace is taken over the ci-electrons, and ft = 
■Ap, The trace is calculated from the eigenvalues A.; (i 
1 ■ - - Y : of the matrix h (first term in H (Eq.(2)). 
The partition function can, therefore, be recast in the 
form 

N 

Z = Y[(J2 er^f^) JJ(1 + e-^M-ri) (4) 

i LUi— 0,1 j 

The corresponding total internal energy is 

E (") = Y f n , + NfEf (5) 

i L J 

The thermodynamic quantities are calculated by the 
averages with the statistical weights, 

AM = J-e-^M ( 6 ) 

of localized f-electron configurations u>, respectively. 

In the simulations, one starts with random configura- 
tion (of /—electrons) at high temperature and then the 
temperature is ramped down slowly until a stable ground 
state is obtained at very low temperature. Then this 
ground state configuration is chosen as the initial config- 
uration for a particular finite temperature. For a fixed 
finite temperature a large number of Monte Carlo steps 
are carried out to make the system reach the equilibrium. 
We have observed that around 1000 Monte Carlo steps 
are sufficient for this. Then the ensemble average of the 
thermodynamic quantity of interest is calculated on next 
(~ 4000) Monte Carlo steps. Now an important question 
arises whether the system has reached equilibrium or not. 
To ensure this we rerun the simulation starting from two 
different initial configurations; firstly from a fully ordered 
configuration of /—electrons and secondly from random 
configuration of /— electrons. After a certain number 
of Monte Carlo steps (around 1000 steps) the internal 
energies for both the cases approach same value. This 
implies that the system has reached equilibrium. A sim- 
ilar approach was earlier employed in the study of phase 
transitions in FKM on a square lattic o 17 ' 29 . 

In order to determine the type of phase transition, we 
have calculated temperature dependence of the charge 
susceptibility \ and the specific heat C v . The % and 
C v are calculated by the fluctuation-dissipation theorem 
(FDT). The \ is related through the FDT to the first or- 
der variance of the density-density correlation function^ 7 -, 
given as, 

X = -^{(gl=i) - (g„=i> 2 ), (7) 

where (...) represent average over different configura- 
tions. The g n is the density-density correlation function 
defined as, 

1 N 

g„ = ^E 2J u(ri)u(r i +T 1 a x + T 2 ay), (8) 

i ri,T2=±n 
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where w(rj) = Ui, a x and a y are the unitary vectors along 
the x and y directions respectively (here we assumed lat- 
tice constant as unity). 

The specific heat (C„) is related through the FDT to 
the internal energy and defined as, 



(9) 



The specific heat and charge susceptibility are calcu- 
lated by different approaches. The specific heat is calcu- 
lated from eigenvalues of the d— electron Hamiltonian. 
The d— electron spectrum changes with temperature. 
Correspondingly total internal energy of d— electrons 
changes with temperature. This change in total inter- 
nal energy of the system with temperature is nothing 
but temperature dependence of the specific heat. It has 
been seen that in a particular temperature range a few 
eigenvalues in the lower part of d— band increase with 
increase in temperature (while other eigenvalues are re- 
maining unaffected). This manifest in the increase of 
d— electron specific heat around that temperature. The 
charge susceptibility is calculated from density-density 
correlation function (<?„). The g n depends on /—electron 
configuration only. 

In addition, we have calculated the thermal average 
(or the ensemble average) of the /—electron occupation 
(wi(T)) on each site (i = 1 to N) and probability energy 
distribution at a temperature T (P(E,T) or simply P(E)). 
The u>i(T) on each site (i) is given as, 



uJi(T) = {u>i) 



£\ e -/3iJ({^}) 



(10) 



where, summation is taken over the different Monte Carlo 
simulation steps at a finite temperature. We have cho- 
sen a random configuration of /—electrons initially. At 
T = an ordered ground state configuration of /— elec- 
trons is obtained. At a finite temperature T, this ordered 
configuration is chosen as initial configuration. The en- 
semble average of the /—electron occupation, cjj(T), is 
calculated for ~ 4000 Monte Carlo steps, after the sys- 
tem reached equilibrium. 

P(E) is defined as the probability of energy E occurring 
at temperature T in the Monte-Carlo simulation^. It is 
equal to the thermal average of the density of states: 



P(E,T) = 



n{E)e 



-E/kT 



Y,E' n{E')e- E '/ kT '' 



(11) 



where n(E) is the density of states for energy E (the 
number of configurations with energy E) and T is the 
temperature of the simulation. 



III. RESULTS AND DISCUSSION 

At finite temperature, only few exact calculation re- 
sults are available on the bipartite (square) lattice with- 
outii and withal correlated hopping. Without correlated 



hopping (usual FKM) and for half-filled case, they ob- 
served that for all finite Coulomb correlation U > the 
ground state configuration is checkerboard type that per- 
sist up to a critical temperature T C (U). At T c the system 
undergoes a phase transition to the homogeneous phase. 
The phase transition is of first order for U < 1 and of 
second order for U > 1 17 ' 31 . With correlated hopping t' 
and U = 0.5, low temperature phase is an ordered phase 
and mainly type of checkerboard, axial striped and seg- 
regated phase. For U = 0.5, and for —1.0 < t' < 0.55 
the phase transition is of the first order and for t' > 0.55, 
phase transition is of the second order. In the regime of t', 
where phase is of the checkerboard type critical temper- 
ature reduces with increasing t' and in the regime where 
phase is of the type of axial striped and segregated, the 
critical temperature enhanced with increase in the t' . 

Perturbative results for large U indicate that to or- 
der 1/U, the FKM can be mapped onto an Ising 
antifcrromagnet (AFM) in a magnetic field, i? e // = 

Y^Uj) W SiS J + + Ef) J2i s i + constant terms, where 
Si = 2u!i — 1, Si = —1,1. The Ising AFM state on a trian- 
gular lattice is frustrated and leads to large degeneracies 
at low temperature. It turns out that this frustration is 
lifted^ in higher order perturbation in 1/U. It is there- 
fore quite intriguing that one would expect the effects 
of frustration to bear on the finite temperature states 
as on-site Coulomb correlation (U) between d— and /— 
electrons and the chemical potential (/i) are varied. 

There is hardly any exact result for the FKM on the 
non-bipartite lattice at finite temperature^. Therefore, 
it is quite interesting to observe the temperature induced 
phase transitions in the FKM on a triangular lattice. We 
have studied the finite temperature induced phase tran- 
sitions as a function of a range of values of U and fillings 
(i) Hf = Nf/N = j and (ii) n/ = | (N4 is constrained to 
N — Nf) of /—electrons using the method outlined in the 
previous section. We look at the effect of correlated hop- 
ping of d-electrons on phase transitions for several values 
of t' e [-1,1] at a fixed U. 



nf 



In Fig. 1(a) we present the variation of specific heat 
(C v ) and charge susceptibility (x) with temperature for 
t'=-1.0 and C/=l. We observe two peak structure in C v 
and X; one at lower and other at higher temperature. Av- 
erage /—electron occupation on each site (o;,-(T)) given 
by Eq.10 at t'=-1.0 and U = 1 for different temperatures 
is shown in Figs.l (b), (c) (d) and (e). The ground state 
configuration of f-electron is regular. This behavior of the 
/—electron configuration persists upto a finite tempera- 
ture ~ 0.7. Above this temperature the /—electron occu- 
pation starts to smear out. At the temperature, known 
as critical temperature T c (^ 1.2), there is a second peak 
in C v and x, which really corresponds to the phase tran- 
sition from regular phase to the homogeneous phase of 
f-electron occupation. 
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(b) T/t= 0.0005 
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(c) T/t= 1.0 



(d)T/t= 1.1 
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FIG. 1: The specific heat C v and charge susceptibility \ ( a ) an d average f-electron occupation u)i(T) on each site ((b) to (e)) 
for different temperatures at U=l, t'=-l.Q and nj = j. Strength of the filled circles are shown in right side of upper panel. 
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FIG. 2: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation L0i(T) on each site ((b) to (e)) 
for different temperatures at U=l, t'—-0.5 and n/ = j. Strength of the filled circles are shown in right side of upper panel. 



In Fig. 2 we have shown the above mentioned quantities 
for correlated hopping f'=-0.5 and U = 1. The ground 
state configuration is again regular. Two peak structure 
is observed in C v and x- We observe that the critical 
temperature T c , decreases with increasing t'. Behavior 
of both quantities C v and x suggest the second order 
phase transition (from low temperature ordered phase to 
high temperature homogeneous phase). 

In Fig. 3 we have shown the C v , x and 0Ji{T) at different 
temperatures for the t'=0 (usual FKM2) and U=l. The 
probability energy distribution P(E) is given in Fig. 4 at 



two different temperatures one at the critical and other 
below the critical temperature. We observe a sharp peak 
at T c ~ 0.6 in C v and x which is directly related to the 
first order transition of the regular phase to the homo- 
geneous phase. These features are clearly visible in the 
variation of uii(T) at different temperatures. We observe 
two peak structure in the P(E) at critical temperature T c . 
If we decrease or increase the value of temperature two 
peak structure disappears and a single peak structure is 
observed in P(E). These features are in accordance with 
the theory of Challa, Landau and Binder—, indicating 
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FIG. 3: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation u)i(T) on each site ((b) to(e)) 
for different temperatures at U=l, t'—Q and n/ = 4. Strength of the filled circles are shown in right side of upper panel. 



0.01 




FIG. 4: Variation of probability energy distribution P(E) for 
different temperatures ((a) to (b)) at U=l, t'=0 and n/ = |. 



X and u>i(T) is shown in Fig. 5 and Fig. 6 for t'=0.5 and 1 
respectively. Two peak structure is observed in C v and \- 
There is second order phase transition from segregated 
phase to homogeneous phase at T c , which corresponds 
to the second peak in the C v and x- I 11 this region the 
critical temperature T c , where phase transition occurs, 
increases with increasing the correlated hopping t' . 

In the usual FKM on bipartite lattice, second order 
phase transition occurs for U > l 17 i 31 . In order to clarify 
the nature of the phase transition in the usual FKM on 
non-bipartite lattice for U > 1, we have calculated tem- 
perature dependence of the C v , x an( i Ui(T) at U = 5 
and have shown in Fig. 7. We observe multi-peak struc- 
ture in C v and x- The ground state configuration, in 
the form of regular phase, persist upto a finite tempera- 
ture. The second order phase transition from regular to 
homogeneous disordered phase occurs at T c = 2.3, and 
corresponds to the last peak in the C v and x- The T c here 
is very large in comparison to that at U = 1 (T c = 0.6). 
Due to large intra-site Coulomb correlation between el- 
and /—electrons, large temperature is required to smear 
the occupation of /—electrons. 

To observe the effect oit' on phase transition at U = 5, 
we have calculated C v , x and uJi(T) at t' G [—1, 1] (shown 
in Fig. 8, for t' = —1.0 only). We observe second order 
phase transition at all t' values. The phase transition 
occurs at a higher value of critical temperature in com- 
parison to that at U = 1. 



B. rif = | case: 

the first order phase transition at T c . With the increas- 
ing value of t' ~ 0.4, the ground state configuration is Three different ground state configurations namely 
the segregated (also metallic) phased. Variation of C v , bounded, axial striped and segregated phases are ob- 
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(d) T/t= 0.6 
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FIG. 5: The specific heat C v and charge susceptibility \ ( a ) an d average f-electron occupation uJi(T) on each site ((b) to (e)) 
for different temperatures at U=l, £'=0.5 and rif — \. Strength of the filled circles are shown in right side of upper panel. 
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FIG. 6: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation Ui(T) on each site ((b) to (e)) 
for different temperatures at U=l, if =1.0 and n/ = -g. Strength of the filled circles are shown in right side of upper panel. 



served^ at U = 1 and t' G [—1,1]. In order to observe 
the stability of these phases at finite temperatures and to 
see the effect of filling of /—electron on T c , we have cal- 
culated C v , x and Wi(T) at {/ = 1 for t' G [—1,1]. Shown 
in Figs. 9, 10, 11 and 12, are variation of C v , \ and u>i(T) 
for U = 1 and t'=—l, 0, 0.5 and 1, respectively. 

The ground state /—electron configuration for t'=-1.0 
to ~ 0.4 and U — 1, is bounded. At finite temperature, 
for all negative values of t', two peak structure is ob- 
served in C v and x ( see Fig9(a) for t'=-1.0). Correspond- 
ing variation of a;, (T) at different temperatures is shown 



in Figs. 9(b), (c), (d) and (e). We observe second order 
phase transition from low temperature bounded phase to 
high temperature homogeneous phase at T c . 

In usual FKM (t'=0), we observe single peak structure 
in C v and \ (Fig. 10(a)). First order phase transition 
is observed at T c =0.45. Axial striped phase is observed 
at about £'=0.5. This phase is not present at i-filling. 
Single peak is observed in C v and x (Fig-11). First or- 
der phase transition from low temperature axial striped 
phase to finite temperature homogeneous phase is ob- 
served about T c =0.275. Phase segregation in the (n f—t') 
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FIG. 7: The specific heat C v and charge susceptibility \ ( a ) an< l average f-electron occupation u)i(T) on each site ((b) to (e)) 
for different temperatures at U=5, t'—Q and rif = j. Strength of the filled circles are shown in right side of upper panel. 
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FIG. 8: The specific heat C v and charge susceptibility \ ( a ) an< l average f-electron occupation uJi(T) on each site ((b) to (e)) 
for different temperatures at U=5, t'=-l.Q and n/ = j. Strength of the filled circles are shown in right side of upper panel. 



plane is observed for t' > 0.6^. Variation of C v , x an( i 
LOi(T) for t'=l and U = 1 at different temperatures are 
shown in Fig. 12. Multi-peak structure is observed in C v 
and x- Phase transition of the second order, low tem- 
perature segregated to high temperature homogeneous 
phase is observed at T c =2.3, which corresponds to the 
last peak in C v and x- The T c is very large in segregated 
phase in comparison to the other phases. 

We have calculated C v , x an d Wi(T) for U = 5 and 
t' G [—1,1] (not shown here) . Multi-peaks are observed in 
C v and x- At all values of t' second order phase transition 



from low temperature ordered phase to high temperature 
disordered phase about T c (much larger than at U = 1), 
corresponding to the last peak in C v and x- is observed. 

Shown in Fig. 13 is the variation of critical temperature 
T c with correlated hopping t' at U = 1 and nf = \ and 
\. We observe, critical temperature decreases in regular, 
striped and bounded phases and increases in segregated 
phase for both fillings. The critical temperature is very 
large in the phase segregation in comparison to regular, 
bounded and striped phases. 

Let us now discuss the numerical results obtained 
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FIG. 9: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation uJi(T) on each site ((b) to (e)) 
for different temperatures at U=l, i'=-1.0 and nj = 5. Strength of the filled circles are shown in right side of upper panel. 
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FIG. 10: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation oj;(T) on each site ((b) to (e)) 
for different temperatures at U=l, t'=0 and n/ = |. Strength of the filled circles are shown in right side of upper panel. 



above in the context of the results obtained on the bi- 
partite lattic e 17 ' 31 at U — 1, different values of if and 
temperature. Important differences in numerical results 
of the spinless FKM on square and triangular lattices 
with correlated hopping are that phase transition occurs 
at very high temperature (around 10 times larger) in tri- 
angular lattice than in square lattice. Phase transition 
is of first order for negative and few positive values of t' 
in square lattice while in triangular lattice it is of second 
order for all negative values of t'. In both type of lattices, 
we observe that second order phase transition occurs at 



large T c than the first order phase transition. In low tem- 
perature regime, the C v and \ decreases slowly for second 
order phase transition and they decrease abruptly for first 
order phase transition with decreasing the temperature. 
In triangular lattice, the critical temperature in regular, 
striped and bounded phases is larger for i-filling in com- 
parison to at ^-filling of f-electrons while nearly equal in 
segregated phase. 

In nut shell, we have studied the phase transitions in 
the spinless FKM on a triangular lattice with correlated 
hopping t'. We have found that the ground state or- 
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FIG. 11: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation u)i(T) on each site ((b) to (e)) 
for different temperatures at U=l, t'=0.5 and n/ = |. Strength of the filled circles are shown in right side of upper panel. 
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FIG. 12: The specific heat C v and charge susceptibility \ (a) and average f-electron occupation cj;(T) on each site ((b) to (e)) 
for different temperatures at U=l, £'=1.0 and n/ = i. Strength of the filled circles are shown in right side of upper panel. 



dered configuration persist up to a finite critical tem- 
perature (T c ). Above critical temperature, homogeneous 
disordered phase is observed. At U = 1 and ^-filling of 
/—electrons, second order phase transition from regular 
to disordered phase is observed for negative values of t' . 
However, for < t' < 0.4 the same phase transition is 
seen to be of first order. Again a second order phase 
transition is observed for t' > 0.4 (from segregated phase 
to disordered). The critical temperature T c , decreases 
in regular phase and increases in segregated phase with 
t'. At i-filling of /—electrons second order phase transi- 



tion is observed for negative values of t' (bounded phase 
to disordered). First order phase transition is observed 
for < t' < 0.55 (bounded and axial striped phase to 
disordered). Second order phase transition is observed 
for t 1 > 0.6 (segregated phase to disordered). The T c 
decreases in bounded and in axial striped phase and in- 
creases in segregated phase. At half filling, segregated 
phase exists at higher values of t' than at one fourth fill- 
ing of /—electron. Correspondingly, T c decreases up to 
a higher value of t' in half filling than one fourth fill- 
ing of /—electrons. The T c , in regular and in bounded 
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FIG. 13: Variation of critical temperature T c with correlated hopping t' for two different fillings of f-electrons at U=l and 
N = 12 x 12. 



phase is larger and similar in segregated phase for one 
fourth filling, in comparison to half filling of /—electrons. 
Second order phase transition is found for all calculated 
values of t' and filling of f-electrons at U = 5. As men- 
tioned earlier there is rarely any numerical calculation 
results available for the spinless FKM on triangular sys- 
tems at finite temperatures. Our results may be helpful 



and motivate to experimental and theoretical study for 
finite temperature properties of the frustrated systems 
of recent interest like transition-metal dichalcogcnides, 
cobaltatcs. Gdl2, NaTi02 and NaV02 etc. 
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